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Notation 


English Letters 
C x , C 2 , C 3 , C 4 , C 5 

^iw> ^2 W 

c 

M 

c s 

C £i * C 62 

D 

1*1 j 1*2 9 1*3 

hi 

J = 0 or 1 

k = ^“(u 2 + v 2 + w 2 ) 

P 

r 

R 

Re = u*R/v 
R d = U 0 D/v 
u 

u * 

U 

U 0 

v 

w 

W 

x 

Y + = (R-r)u*/v 


model constants, specified in Table 1. 

model constants introduced by wall 
correction, specified in Table 1. 

model constant associated with v t , 
specified in Table 1. 

model constant associated with diffusion 
model, specified in Table 1. 

model constants associated with 
r-equation, specified in Table 1. 

pipe diameter 

damping functions, specified in Table 1. 
grid spacing 

index denoting two-dimensional or 
axi symmetric flows . 

turbulent kinetic energy 

fluctuating pressure 
radial coordinate 
pipe radius 

turbulent Reynolds number 
pipe Reynolds number 

fluctuating velocity along x-direction 
friction velocity 

mean velocity along x-direction 
mean velocity at pipe center 
fluctuating velocity along r-direction 
fluctuating velocity along ^-direction 
mean velocity along 0-direction 
axial coordinate 
normalized r-coordinate 


v 



Greek Letters 


i 


a = 1/k 


P 


a l 9 P\ i 


C 

B 

K 

V 



P 

O’ 

€ 


Q 


log-law slope 

constant in law- of- the- wall 

model constants associated with 
redistribution model, specified in 
Table 1. 

dissipation rate of k 
circumferential coordinate 
von Karman constant 
fluid kinematic viscosity 
turbulent diffusivity 
fluid density 

model constant associated with 
e-equation, specified in Table 1, 

normalized dependent variables 

angular velocity 
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Abstract 


A full Reynolds-stress closure that is capable of describing 
the flow all the way to the wall has been formulated for 
turbulent flow through circular pipes. Since viscosity does not 
appear explicitly in the pressure redistribution terms, 
conventional high-Reynolds-number models for these terras are 
found to be applicable. However, the models for turbulent 
diffusion and viscous dissipation have to be modified to account 
for viscous diffusion near a wall. Thus modified, viscous 
dissipation in the flow is no longer isotropic as postulated by 
Kolmogorov for high-Reynolds-number turbulence. Two 
redistribution and two diffusion models are investigated for 
their effects on the model calculations. Wall correction to 
pressure redistribution modelling is also examined. Diffusion 
effects on calculated turbulent properties are further 
investigated by simplifying the transport equations to algebraic 
equations for the Reynolds stresses. Two approximations are 
explored. These are the equilibrium and non-equilibrium 
turbulence assumptions. Finally, the two-equation closure is 
also used to calculate the flow in question and the results 
compared with all the Other model calculations. 

Fully-developed pipe flows at two moderate Reynolds numbers 
are used to validate these model calculations. They are chosen 
because detailed turbulence measurements near the wall are 
available. The calculations show that all closure models give 
good agreement with measurements of mean velocity, shear stress, 
turbulent kinetic energy and dissipation rate near a wall. 


vii 



However, the slope of the logarithmic law-of-the-wal 1 recovered 
from these calculations varies from one closure model to another. 
Some closure model predicts the correct behavior of the log-law 
constant as a function of Reynolds number, while others provide 
the wrong trend. Wall correction is found to have little effect 
on the model calculations. Mean-strain effects on redistribution 
modelling are found to give rise to an adverse influence on the 
calculated log-law, in the case of non-equilibrium algebraic 
stress closure. All closure models examined fail to predict the 
steep rise of turbulence intensities near a wall correctly. 

Also, they fail to reproduce the isotropic behaviour of the 
normal stresses at the pipe center. Overall, the best model 
prediction is given by the full Reynolds-s tress closure 
incorporating a non-isotropic gradient diffusion model and the 
Launder et al. (1975) model for pressure redistribution. 
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I n t r oduc t i on 


1 . 

1 . 1 Mo tivation 

The flow inside a practical combustor, such as gas turbine 
or solid fuel ramjet combustor, is usually very complicated, 
especially at the combustor wall. Besides geometry effects, 
cooling and dilution air through the combustor wall, fuel 
sublimation from the solid fuel grain and rotation of the 
combustor further complicate the wall boundary- 1 ayer flow. 
Conventional modelling of combustor flow is to assume high- 
Revnolds number turbulence, since the flow Reynolds number is 
normally very high inside the combustor. The wall boundary 
conditions are satisfied indirectly by specifying some empirical 
functions, such as the logarithmic law-of-the-wal 1 for mean 
velocity and equilibrium turbulence for the turbulence field, to 
link the conditions at the wall to the first calculation point 
away from the wall. Consequently, the effects of viscosity and 
conductivity near the combustor wall cannot be resolved correctly 
and hence their influence on the flow outside the near wall 
region cannot be assessed. Therefore, this presents a difficult 
problem for the calculation of flow inside practical combustor 
because the complicated boundary conditions render the simple? 
law-of-the-wal 1 and equilibrium turbulence assumptions near the 
wall invalid. It is clear that conventional hi gh~Reyno 1 ds-number 
turbulent closure models need to be modified to account for 
viscosity and conductivity effects ne?ar a wall before they can be 
applied with confidence to calculate the flow inside practical 
combustors . 
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Some work along this direction has been carried out by 
previous researchers (e.g. Jones and Launder 1972b; Han jail e and 
Launder 1976; Chien 1980). These studies will be reviewed in the 
next section and further improvements will be identified. 
Therefore, based on this review, the objectives of the present 
study are) formulated in Section 1.3. 

1.2 Brief Discuss ion of Previous Work 

Reynoldses tress closure + of turbulence applied to turbulent 
flow calculations was first examined bv Hanjalic and Ltiunder 
(1972). In their model, certain assumptions concerning the 
structure parameters, u i u i /k, (i not summed), where u i is the i th 
c omp o n o n t of the fluctuating velocity and 2k - u . u x ( s umm a t i o n 
over i) is the turbulent kinetic energy, were invoked to simplify 
the four Reyno 1 ds-s t ress transport equations for two-dimensional 
thin shear layers to two equations for- turbulent shear stress and 
k. These were then solved with the mean flow equations and an 
equation governing the transport of <•: , the dissipation rate of k. 
The closure was arrived at by assuming the flow Reynolds number 
to be very large and that u i u i /k - constant throughout the shear- 
layer with the constants given by plane shear flow measurements. 

In view of these approximations, the boundary conditions cannot 
be applied at the wail. Rather, they were applied near the wall . 
Tn particular, the mean flow velocity was matched to the 
logarithmic law-of-the-wal 1 , the gradient of k was set equal to 

+ The term " Rev n o I ds-s t ress closure” is used to done to closure 
schemes that solves the full set of Revno 1 ds-s tress transport 
equations as well as models that solve the Reynolds shear stress 
and k equations alone ( e . g . Hanjalic and Launder 1972, 1976). 
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zero, the shear stress was determined from the mean momentum 
equation with the convection terms neglected and r was set equal 
to the turbulence generation rate. The model gave good 
comparison with measurements away from the wall for a wide 
variety of thin shear layers. However, detailed flow modelling 
near a wall remained unattainable, just as in the case of the 
mixing- length model where the same equilibrium turbulence 
arguments were used to determine the behaviour of the mixing 
length near a wall. 

In view of the initial success of the Revno 1 ds-s t ress 
closure model, later researchers (e.g. Launder et al. 1972; 

Mellor and Herring 1973; Mellor and Yamada 1974; Irwin and Smith 
1975; Launder et al. 1976; Gibson and Rodi 1981) relaxed the 
assumption, u i u i /k constant throughout the shear layer, and 
solved the full set of Reynolds-stress transport equations. 
However, the large-Revnolds-number assumptions were retained in 
the modelling of the turbulent diffusion, redistribution and 
energy dissipation rate terms in the Reyno 1 ds-s t ress transport 
equations. Consequently, the boundary-layer flow very near the 
wall had to be handled in the same manner as that proposed by 
Hanjalic and Launder (1972). The boundary conditions for u i u i 
near a wall, however, required special attention. In general, 
either a slip condition for u i u i was imposed (Irwin and Smith 
1975) or the Neuman boundary conditions were specified (Mellor 
and Yamada 1974). The amount of slip specified for u i u i depended 
to a great extent on the type of flow considered. As a result, 
the closure model was problem dependent and, in spite of the 
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improvement, it still could not provide an accurate description 

of the flow very near the wall. This, in turn, means that the 

/ 

Revnolds-st ress closure cannot be used to estimate the Reynolds 
number effects on turbulent flows, because in the immediate 
vicinity of a wall viscous effects have to be important (Mellor 
and Herring 1973). 

The logarithmic law~of- the-wal 1 assumption was generally 
applicable for a wide class of simple turbulent flows. However, 
it failed to provide a reasonably accurate estimate of the near 
wall mean velocity in separating flows, in r e 1 am in ar i z ing flows 
and in complex turbulent flows (e.g. Stratford 1959; Jones and 
Launder 1972a; Bissonnette and Mellor 1974). In order to remedy 
this unsatisfactory boundary condition, Jones and Launder (1972b) 
proposed to modify the two-equation model of turbulence, i.e. the 
k-c model, for Reyno 1 ds-nurab er effects. The shear stress, -u 1 u 2 , 
was then calculated by assuming - u l u ? = v t (0U 1 /dx 2 ), where 
v t = C k 2 / e , Uj was the local mean flow and x 2 was the 
coordinate normal to the wall. With suitable modification to c 
near a wall to account for viscosity effects, they found that the 
modified k-c equations can be integrated with the boundary 
conditions, U 1 = 0, k-0 and e-0, applied at the wall. This 

allowed the near wall flow to be calculated directly from the 
governing equations, and good agreement with r e 1 am i n ar i z i ng flow 
measurements (Jones and Launder 1972a) was obtained. 

Furthermore, the logarithmic behaviour of the near wall flow was 
recovered when the flow Reynolds number was sufficiently large. 
Later, Hanjalic and Launder (1976) applied the arguments of Jones 
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and Launder (1972b) to modify their (Hanjalic and Launder 1972) 
Reynolds-s t ress closure model to account for viscosity effects 
near a wall. In their new closure model, they relaxed the 
assumption Uj u t /k - constant. Instead, they assumed 
(3/4)(u 2 + u| ) - k and u| = 4(u 1 u 2 ) 2 /k based on the pipe and 

channel flow measurements of Laufer (1954) and Eckelman (1970). 
This way, the near wall behaviour of u 2 /k, i.e. u?/k 0 as 
x ? -* 0, was satisfied. Their calculated mean velocity and shear 
stress results were in excellent agreement with the channel flow 
measurements of Patel and Head (1968) and Eckelman (1970) and the 
relaminariz.ing flow data of Jones and Launder (1972a). In spite 
of these successes, the near wall behaviour of u, u. /k cannot be 
calculated. To do this, one needs to resort to a full Reynolds- 
stress closure model where all the transport equations for iru. 
are solved rather than the equations for k and \i l u 2 alone. 

1.3 Present Objectives 

The primary objective of this study is to formulate a full 
Reyn o 1 ds-s t res s closure model so that the calculations can be 
carried all the way to the wall and satisfy the boundary 
conditions at the wall for U, , ^ and c. Validation of the 

model is carried out by comparing the calculated results with 
fully-developed turbulent pipe flow data at two moderate Reynolds 
number. By selecting fully-developed turbulent pipe flows, the 
complexity involved in solving the transport equations can he 
greatly reduced because the governing equations simplify to 
second-order, non-linear ordinary differential equations. This, 
in turn, allows the various modelling assumptions to be assessed 
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easily. A secondary objective, is to investigate the effects of 
redistribution and diffusion models on the modelled flow. To 
this end, the models of Rotta (1951) and Launder el a 1 . (1975) 

for pressure redistributions are examined in detail together with 
a non- isotropic and an isotropic gradient diffusion model for 
turbulent diffusion. The effects of diffusion modelling are 
further investigated by simplifying the transport terms in the 
Reynolds-stress equations according to the suggestion of Rodi 
(1976). This results in a set of algebraic equations for the 
Reynolds stresses and can be solved with the 1 ow~ Reyn o I ds~ numb er 
form of the k-c equations (Jones and Launder 1972b; Chien 1980). 
The final objective, then, is to compare all the above model 
calculations with the results obtained from the basic two- 
equation model and to identify a model thaf performs the best :i n 
pipe flow calculations. 
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2 . Low- Reyn olds -Numbe r Closu res 

For an incompressible flow, the transport equations for the 
Reynolds stresses u f u^ can be concisely expressed in Cartesian 
tensor as 


i” + 5 jk u i ) 


t d u, u . 

4 - y — — > 

a *k J 


au ; , au, 

u i u lf-S — + Uj U V T — L 

J R dx k ’ K dx k 


+ 

/» ax j 


-- 2 v 


du i <3u • 

ax k ax k 


(i) 


Here, lower and upper case ,u’s denote fluctuating and time- 
averaged velocity components, respec t i ve ly , and overbars imply 
the usual time averaging of the correlations in question. 

The terms in (1), from left to right, in general, can be 
interpreted as the convection, diffusion, production, 
redistribution and viscous dissipation of n i u^ , respectively. Of 
these five groups of terms, the convection and production terms 
are exact and do not need modelling. In the past, only high- 
Rey n o 1 ds-n umbe r models have been proposed for the diffusion, 
redistribution and dissipation terms. Consequently, the 
resultant closure model is not valid for flows near a wall 


(Mellor and Herring 1973). Although some advances toward this 
direction have been made by Hanjalic and Launder (1976), a full 
Reyno 1 ds-s t r ess closure for 1 ow-Reynol ds-number turbulence is 
still not available. In the next section, an attempt will be 
made to close (1) so that the resultant transport equations are 
valid for low as well as hi gh-Reynol ds-number flows. At least 
two different models are proposed for each of the three terms 
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that required modelling. The relative merits of these models 
will be investigated. Diffusion modelling is further examined by 
greatly simplifying the transport equations into algebraic 
equations for the Reynolds stresses. Two approximations will be 
investigated; one is the equilibrium turbulence assumption and 
another is Rodi’s (1976) approximation. A discussion of these 
algebraic stress closures is given in Section 2.2. Finally, the 
two-equation closure model of Chien (1980) for 1 ow-Revno 1 ds - 
number turbulence is included in Section 2.3 for the sake of 
completeness . 

2 . 1 Full Revnol ds-S tress Closures 

In order to model (1) for 1 ow-Revno lds-number turbulence, 
appropriate models for the redistribution, diffusion and 
dissipation terms have to be formulated. The subsequent sections 
provide a first attempt for this endeavour. 


2.1.1 Redistribution Models 
Since the term 


P 

P 


(5Hj- + 

ox . 


ax, ; 


has a zero trace for an incompressible flow, it acts to diminish 
the difference between the normal-stress components (Hinze 1959). 
Therefore, it neither produces nor destroys turbulence energy. 
Furthermore, p satisfies the equation 


1 d 2 P _ _ f d 2 (Uj Uj ~ Uj Uj ) + q 3u, 1 

p 0X S 0Xj 0Xj 0X. 0X • 0 X . -I * 


( 2 ) 
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obtained by taking the divergence of the equation for . Since 
v does not appear explicitly in (2), this suggests that, to first 
order, any h i gh Revn ol ds -number model for this term can be 
adopted for the present study. Specifically, the redistribution 
model proposed by Launder et al. (1975) is adopted. This can be 
written as 


p K dx i dx i } 


where 


C ik (Ui U -* ~ 3*i i k) 

“ “ /»i(D u 

- iCj k S ■ . , 


|^ P ) 




[ u j u 


3U ; 


au ; 


k dxt + U ‘ Uk 0x L 


]. 


D iJ 


r au„ 

= - L u i u k ^ + u j u 


k 0X ; 


au k 
‘J “ k ax,. 


]. 


s i. 


p = - 


«Hj. + 22j. 

6Xj dXj ’ 

dU. 

u< U b L- 


l i 


dx. 


(3) 


and C 1 , a 1 , p l , and , are model constants. According to 

Launder et al. (1975), a x , /> T , and *\ are not independent 
constants. Rather, they are related to one constant C 2 . These 
and other model constants are listed in Table 1 for reference. 

Launder et al.’s model includes both the symmetric and 
antisymmetric mean-strain effects on redistribution modelling. 
However, their influence on near wall flow calculations has not 
been clearly demonstrated. In order to evaluate the effects of 
mean-strain modelling on the flow near a wall, Rotta’s (1951) 
simple return-to-isotropy model for the redistribution terms will 
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also be examined in the present investigation. The model is 
given by (3) by setting = & 1 = - 0 . 

2.1.2 Wall Correction to Redistribution Modelling 

Since the presence of a rigid wall affects the pressure 
field, thus impeding the transfer of turbulence energy from the 
streamwise direction to that normal to the wall, Launder et el. 
(1975) propose a wall correction to the pressure redistribution 
model to account for this wall effect. The correction is 
designed specifically to model the decrease of turbulence energy 
transfer to the normal direction. Since then, the wall 
correction has been used by Irwin and Smith (1975) to model 
curved shear flow and by Gibson and Launder (1978) to model 
atmospheric boundary layers. However, in these calculations, the 
near wall flow is not resolved directly. Therefore, the value of 
wall correction in pressure redistribution modelling has not been 
clearly demonstrated. The present approach allows the near wall 
flow to be calculated directly and, thus, provides a good 
opportunity to assess t h e relative merits of the w a 1 I corr e c t i o n . 

In view of this, the wall correction proposed by Launder 
et al . (1975) will also be investigated. When this correction 

term is included in the pressure redistribution modelling, the 
comp lete model becomes 
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d , du ■ du ■ 

- (4— > ~ + 


P 'dx,, 


dx 


- c i £( u > u j - 3 5 i j k > 


a 1 ( ^ i j “ -3 ^ i j ^ ) P 1 ( ® i j 3^ij^ 


- k S t , 


[C lw rCUiUj - jkjP) 


+ C 2W (P U - D j j ) ] , 


(4) 


where x 2 is measured normal to the wall and the model constants 
C }w and C 2w are specified in Table 1. 

2.1.3 Diffusion Models 

The other terms in (1) that need modelling are the diffusion 
and viscous dissipation terms. Since these terms involve v 
explicitly, their hi gh-Reyno lds-number models (Launder et al. 
1975; Kolmogovor 1941) have to be modified to account for 
viscosity effects. This can be easily carried out for the 
diffusion model by including the term v(3u. u^ /6x k ) in the final 
diffusion model for low-Revnolds-number flows just as Hanjalic 
and Launder (1976) have done. If pressure diffusion is 
neglected, as suggested by Hanjalic and Launder (1976), the model 
for the diffusion term becomes 


[■ u i u j u k + v £ r Ul - + 7(^k u i + 5 ik u j) ] 


dXjj L i j k c)x k P 


r 3u . u . k , (Ju. u . , — 

- + c -7 < ‘."itV u i u >-ai7" 


•du. u. 


■3u. u ; 


dx 


du.u^. -i 

- u * u > _ ax7^ ) J 


(5) 
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This model is based on the gradient diffusion assumption and 
is tensorially correct. Also, it gives a non-isotropic 
diffusivity. On the other hand, past researchers (e.g. Mellor 
and Yamada 1974) have found that a much simpler isotropic 
diffusivity model for high-Reynolds-number turbulence works 
equally well as the non-isotropic one, even though the model is 
tensorially inconsistent. Whether this also holds true for low- 
Reynolds-number turbulence will again be examined in the present 
investigation. To this end, the simpler isotropic gradient 
diffusion model given by 


a 

0X k 


[- u i u j u k 


3U; U 4 

V — — 1 J- 

dx,„ 


:(*» 


jk u i 


+ 5 ik u j)J 



where 


v 


t 


♦ *.> ■ 


is defined as 


c — f, 

fl c * 


( 6 ) 


(7) 


and f 2 is a damping function specified in Table 1, is also 
investigated. It should be pointed out that inherent in this 
model is the assumption that pressure diffusion is not important 
and can be neglected. A comparison of these two models will, 
therefore, provide a clear indication on which gradient diffusion 
model is more appropriate for 1 ow-Reyno 1 ds-numb er turbulence 
closure. 

2.1.4 D i s s i p a t i on Models 

For high-Reynolds-number turbulence, Kolmogorov (1941) 
assumed viscous dissipation to be isotropic. However, near a 
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wall this assumption is no longer valid because the turbulent 
flow Reynolds number in this region is not large. To see how 
Kolmogorov’s dissipation model should be modified to account for 
Reynolds-number effects, the Reynolds-s tress equations (1) are 
examined for near wall behaviour with the proposed diffusion 
models (5) and (6). 

When (5) is used in conjunction with the Revnolds-st ress 
closure of Hanjalic and Launder (1976), good results for channel 
(or two-dimensional) flows are obtained. However, when it is 
used in conjunction with the full Reynolds-s tress closure for 
axisymmetric flows, the following difficulty would appear at the 
wall. The difficulty is associated with the term 
(6/dx k ) (v0u. Uj /0x k ) . If cylindrical coordinates (x, r, 8) are 
used to expand this term for fully-developed pipe flows, then at 
the wall, the leading term becomes 


1 a 

r J 0r 


(r- 


0U. U. x 

V— ^ *-) 

0r 



0 U . U . v 
0r ; 


+ 


J v 0U . u. 
r 0r 


( 8 ) 


where J = 0 or 1 for two-dimensional or axisymmetric flows, 

respectively. At the wall, u. u.. = a i ^.y N + b.,.y N + 1 + (Mellor 

and Herring 1973), where N > 2, ^ , b.^, are constants to be 

determined, and y = R - r. Substituting this expansion into (8) 
gives 

_1_| ir J v *±Jh-) = (N-l ) Naj j vy N “ 2 

+ [N (N+l )b J j V - + O(yH) . (9) 

Therefore, for N=2, all the terms on the right hand side of (9), 
except the first term, vanish at y = 0, and this is true for all 
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near wall flows, be it two-dimensional or axisymmetr ic . This 
means that molecular diffusion is finite at the wall. Since the 
term ( d/dx k ) ( vdu t u^ /dx k ) does not need modelling and v does not 
appear in other diffusion terms; additional terms are required in 
the modelling of the dissipation function, 2 v ( du A /dx k ) ( du^ /dx k ) , 
in order to balance the finite molecular diffusion at the wall. 
The foregoing arguments, therefore, suggests that the Kolmogorov 
(1941) high-Reynolds-number model for the dissipation function 
should be modified to give, 


2 v 


du . du ■ 

ax k ax k 


3 * ijC (* 2^)2 


for low-Reynolds-number flows. 

The dissipation model (10) is not i 
Reynolds-number is finite. However, at 
number (10) approaches Kolmogorov’s mode 
correct limiting behaviour is recovered. 


sotropic when the 
very large Reynolds 
1 asymptotically and 


( 10 ) 


the 


2 . 2 Algebraic Stress (ASM) Closures 

Less sophisticated closure models for high-Reynolds-number 
turbulence have also been put forward by various researchers. 
Specific assumptions are put forward to simplify the Reynolds- 
stress equations (1) so that the equations are reduced to 
algebraic equations for the Reynolds stresses. The equilibrium 
turbulence assumption is used by So (1975, 1977) and So and 

Mellor (1978) to calculate curved shear flows, rotating and 
swirling flows. On the other hand, Rodi (1976) proposes to 
approximate the transport (convective and diffusive) terms in (1) 
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by the stress ratios ^u^/k. Since then, this non-equilibrium 
turbulence scheme has been used by Gibson (1978), Gibson and 
Launder (1978) and Leschziner and Rodi (1981) to calculate a wide 
variety of turbulent shear flows. These closure models will also 
be extended to low-Reynolds-number turbulence and their 
performance compared with the full-Reynolds-s tress models to 
further identify the effects of diffusion modelling. 

In the course of modelling (1) taking the algebraic stress 
closure approach, the solution of two more transport equations 
are required (Gibson 1978). Normally, the k and c equations are 
solved in addition to the mean flow equations. Therefore, the 
modifications of these closure models for low-Reynolds-number 
turbulence can be achieved via two different approaches. One is 
to modify the basic equations for k and c and another is to 
modify both the k-c equations and the models for pressure 
redistribution and dissipation. In Section 2.1, it has been 
shown that the high-Reynol ds-number form of the pressure 
redistribution models is also applicable to low-Reynolds-number 
flows. However, the dissipation model has to be modified for 
low-Reynolds-number turbulence. In algebraic stress closures, 
turbulent diffusion is either neglected or approximated by the 
stress ratios via Rodi * s approximation. If the dissipation model 
is modified to account for viscosity effects, then it can be 
shown that, under the assumption of equilibrium ASM, the 
resultant algebraic stress equations are not balanced at the 
wall. On the other hand, under Rodi’s approximation, the 
additional dissipation term would cancel out with the extra 
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dissipation term introduced by Chien (1980) in the k equation. 
Consequently, the approach taken here is to modify the k —c 
equations alone. The modified forms of the k -c equations are 
given in Section 2.3, while the simplified forms of the 
Reynolds-stress equations needed for the equilibrium and non- 
equilibrium algebraic stress closures are specified in the 
following two sections. 


2.2.1 Equilibrium Assumption 

If the equilibrium turbulence assumption is invoked, 
production of turbulence energy is equal to its dissipation rate. 
The Reynolds-stress equations (1) simplify to 


, dU. -| 

* u i u «c5uH 


+ Ef^SLu + 

p ox. 0X 5 


6u i <3u ■ 

dx k 0^ 


(ID 


The hi gh-Reynol ds-numb er models for the pressure redistribution 
and dissipation terms are used to close (11). These are given by 
(3) and the Kolmogorov isotropic dissipation mode], or 


du. <3u . 2 . 

t 1— 1- ~ A , . £ 

dx k dx- K 3 


Again, Rotta’s re turn- t o-isotropy model is obtained from (3) by 
setting a - p Y = i x = 0. Both Rotta’s and Launder et al.’s 

models will be investigated and the results compared with the 
full Reynolds-stress models. 
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2.2.2 Non-Equilibrium As sump t ion 

In this case, the transport terms in (1) are approximated by 
Rodi’s (1976) assumption and (1) is reduced to 


(P - 


~ C u j u n; 


+ U; U, 


+ 


P, du. + QUj 
P v ax k ax k 


2 v 


0u i du . 

dx k ax k 


(13) 


If (3) and (12) are used to model the last two terms in (13), 
algebraic equations for u A u^ can be derived from (13). It is now 
clear that if (10) is used to model the dissipation term, c in 
(13) will have to be replaced by c + 2vk/x.| (see eq. (14)). 

The net result is again equivalent to (13) with the dissipation 
function given by (12). 


2 . 3 Two-Equation Closure 

Equation (1) can now be expressed in terms of u i u^ , U. , k, e 
and their gradients. Transport equations are, therefore, 
available for all unknowns except e . Two options are available 
for c. Either the equation proposed by Jones and Launder (1972b) 
or the equation modified by Chien (1980) for 1 ow-Reyno 1 ds-number 
turbulence can be used. Since the k equation obtained by 
contracting (1) using the proposed models is similar to that 
proposed by Chien (1980), it is decided that Chien’s c-equation 
would be more appropriate for the present full Reynolds — stress 
closure. Also, (12) and (13) can be expressed in terms of u. , 
U. , k, c and the gradients of U. . To complete the definition of 
k and c in the flow field, two equations governing the transport 
of k and r are required. Since the ^-equation modified by Chien 
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(1980) is used in the full Reynolds-s tress closure it would be 
appropriate to also use the k-equation modified by Chien for low- 
Reynolds-number turbulence rather than by contracting (1) with 
models given by (3), (5) and (10). These two equations are: 


Dk d r, x dk -] „ 

5t = 557 L (v + + p - 


C - 


2vk 


(14) 


Dc 

Dt 


dx h 


[( 


v + 


dc 


dx. 


] + C 6ik P C 62 f ii, 


2 vc ~C 4 x 2 u* /v 


x l 


(15) 


where x 2 is measured normal from the wall and the model constants 
C f i , C c2 , c 4> cr £ and damping function f x are specified in 
Table 1. The diffusivity v t is again taken to be given by (7). 
Therefore, the two-equation closure is given by solving (14) and 
(15) together with the mean flow equations. It should be pointed 
out that, in the full Reynolds-s tress and algebraic stress 
closures, an assumption for the turbulent stresses, u i , is not 
required. For the two-equation closure, a gradient transport 


model for u t Uj is invoked, such that 



3. Ful ly- D evelope d Pip e Flows a t Mo derate Re ynolds Num bers 

Fully-developed turbulent pipe flows at moderate Reynolds 
numbers are used to validate the full Reynolds-s tress and 
algebraic stress models. If models (4), (5) and (10) are 
substituted into (1) and the component equations are written in 
cylindrical coordinates (x, v, 8) with mean and fluctuating 
velocities given by (U, 0, 0) and (u, v, w) , respectively, then 

the governing equations including the simplified mean flow 
(Laufer 1954) and c equations are: 
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dU 

dr 
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u|r 
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( 22 ) 


In writing down these equations, a damping factor f 3 specified in 
Table 1 and suggested by Hanjalic and Launder (1976) is included 
in the stress production term of (22). This damping factor is 
found to be necessary in the course of solving (17)-(22), because 
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without it shear production near a wall is found to be over- 
predicted. Similarly, turbulent di f f us i v i t ies near a wall are 
found to be way over-predicted if C s is taken to be a true 
constant. In all subsequent calculations, C p is damped by a 
factor f 2 similar to that specified in (7). These and other 
model constants are specified in Table 1 and a re consistent with 
those recommended by various researchers (Rotta 1951; Launder 
et al. 1975; Chien 1980). 

Boundary conditions are specified at the wall and at the 
symmetry plane. These are given by: 

U~ e = u 2 -v 2 -w 2 - uv =0 at r = R (wall), (23) 


de _ du 2 _ dv 2 _ dw 2 
dr dr dr dr 

uv = 0 



at r = 0 , 


(24) 


Only one boundary condition is specified for U, since the 
symmetry condition for U has been utilized to evaluate the 
integration constant when the mean momentum equation is 
integrated to give (17). 

Similarly, component equations for (1) using (4) and (6) for 
closure can also be written down. However, for brevity’s sake, 
they are not included here. As before, it is found necessary to 
damp the shear production term in the equation for uv by the 
factor f . Rotta’s model is given by the same set of equations 
with a 1 = - Y 1 - C lw = C 2w = 0. If no wall correction is 

required, C lw and C 2w should be set equal to zero. 
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For algebraic stress and two-equation closures, the 
transport equations to be solved are (17), (18) and a 

corresponding k-equation which can be written as 

2 vk 


1 d r, , x dkq — ,dU. 

F 37 » + -t)^ - “''<dF> - 


(R-r)* 


= 0 


In the case of algebraic stress closures, the stresses u t u^ 
provided by (11) or (13) which can be written out in their 
component forms. For two-equation closure, -uv is provided 
(16) which reduces to 


- uv = v* 


K dr' 


The boundary conditions for k are: 


k = 

dk 

dr 


0 

0 


at r = R, 

at r = 0 . 


( 25 ) 

are 

by 


(26) 


(27) 
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4. 


Method of Solution 


The set of equations (17)-(22) with boundary conditions (23) 
and (24) are solved numerically by the Newton iteration scheme 
(Na 1979). First, normalization of the dependent variables U by 
u* , Uj u. by u 2 and c by u^/R are carried out, while the 
independent variable r is made non-dimensional by u*/v, so that 
Y* = u*(R-r)/v is the new dimensionless coordinate. Therefore, 
the integration from the pipe wall to the centerline is now 
carried out from Y + = 0 to Y + = Re, where Re = u*R/v is the 
turbulent Reynolds number to be specified. Re is related to the 
pipe flow Reynolds number R D - U Q D/v by Re = (u*/2U Q )R D . 
Therefore, once R D and the pressure drop along the pipe are 
known, Re can be determined and it becomes the only input 
parameter to the problem. 

Next, the six first- and second-order ordinary differential 
equations are written into eleven first-order equations by 
defining new variables for d ( u i u^ / u| ) /dY + and d ( c R/u 2 ) /dY + . 

If these eleven variables are denoted by , i = 1, .... 11, 

such that iTj = U/u* , ^ 2 - c R/u 2 , ^ 3 = u 2 /u 2 , = v 2 /u 2 , 

= w 2 /u 2 , = nv/u 2 , y 7 = d^ 2 / dY + , = d* 3 /dY + , = d^ 4 /dY + , 

^ 1Q = d^ 5 /dY + , a = d^ 6 /dY + , then their derivatives can be 
approx i mat ed by cen t ered-d i f f er ence gradients and averages 
centered at the midpoints of the grid, defined by 
Y+ r. 0; Y+. = Y + .. x + h. , j = 1, 2, . . . , M; Y + M = Re; (28) 


d »i - (±±jU z it 

dY+ ~ h j 


(29) 
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The resultant equations are nonlinear algebraic equations. 
Therefore, they have to be linearized before the; factorization 
scheme can be used. Newton’s linearization scheme is used so 
that the Newton (k+1) iterates for (^ t ). can be written as 

(*i)/ k + 1 > - (*, )/ k > + SiVt )j< k > . (30) 

These are then substituted into the eleven first-order 
differential equations. If quadratic and higher-order terms in 
5 ( k > are neglected, the resultant linear algebraic equations 
can be put into v e c t o r ~ matrix form as 


[A] [ 6 ] - [*] , 


(31) 


where [A] is the coefficient matrix of order MxM and its elements 

a r e in a t r i c e s of order 11x11. The matrix [ A ] is of the 

tri diagonal form, while the matrix [ ¥ ] is a column matrix of 


( 1 *. ) and (**. 


j - i 


The boundary conditions are: 


<H>i ) 0 


d (</;■-. ) 

i ) M 


= 0 


, 6 ; 


(32) 


7, 


10 


0 n c e t h e e q u a i i o n s a r e p u 1 into the form { 3 i ) and (32), the v c a n 
be solved in t e rat i v e 1 y us in g a n y mat r i. x L n version technique. 
Iteration i s c a r r ied o u t un t i 1 [ 6 ] me e t s c e r t a :i n a ecu r a o y 

criterion. For the present study, the accuracy criterion is 


chosen as 


) 4 


' - $ ( iP ■ 

d ( y. . ) _ 4 ' 


< I: - : 


< 10 - 4 


(33) 
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A non-uniform grid is used to carry out the calculations. 
Typically, five grid points are specified between Y + = 0 and 
Y + = 5. This is followed by 15 grid points between Y + “ 5 and 
Y + = 65. The rest of the region 65 < Y + < Re is then divided 
into 30-50 grid points depending on the problem considered. In 
general, this system of grid spacing is sufficient to give 
convergent solution after ^1500 iterations. 

The equations for the other closure models can be similarly 
solved. However, for brevity’s sake, they are not outlined in 
detail here. Anyway, once a solution scheme has been developed 
for the more complex set of equations (17)-(22), the same scheme 
can be easily adapted to solve a set of simpler equations. 


25 



5. 


Presentation of Results 


Fully-developed pipe flows at two different Reynolds numbers 
are selected for comparison with the model calculations. These 
are the detailed measurements of Laufer (1954) at R D = 50,000, 
and of Schildknecht et al. (1979) at R D = 21,750. They are 
chosen for their careful measurements of the turbulence field 
near the wall. Consequently, they would provide accurate data 
for the evaluation of the full Reynolds-st ress model for low- 
Reynolds-number turbulence. The input parameter for these two 
calculations is Re = 1052 and 489, respectively, for Laufer ’s and 
Schildknecht et al.’s experiments. Calculations are carried out 
to compare the performance of different redistribution and 
diffusion models. The redistribution models examined are Rotta’s 
(1951) return-t o-isot ropy model and Launder et al.’s models 
including mean-strain and wall effects. As for diffusion, two 
models are investigated; one is Launder et al.’s model given in 
(5) and another is an isotropic model given in (6). The effects 
of diffusion modelling are further examined by considering less 
sophisticated closure models such as algebraic stress closures. 
Again, redistribution modelling effects are studied by comparing 
the Rotta and Launder et al . models. Finally, these mode] 
calculations are compared with the simple two-equation closure 
results. 

In view of the large number of calculated results and their 
close similarities, the comparisons with data are presented 
separately rather than together in one figure for each flow 
property. Although this involves many more figures to be 
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presented, it will give a clearer comparison of each closure 
model with measurements. Altogether eleven model calculations 
are made for each experiment. These are organized for 
presentation in the following manner. With the exception of the 
k -e model calculation, each set of comparison consists of eight 
figures where the model calculations of both experiments are 
presented. These eight figures, numbered a through h, show the 
comparisons of U/u* in semi-log plot, U/U Q , uv/u|, k/u^ , 
cR/u^, u. u A /u^ , u“~Tr. /k and uv/k versus l-2r/D in the sequence 
given. The comparisons of the near wall behavior are shown as 
insets in each figure, and the coordinate used is the normalized 
wall coordinate Y + . Only five figures are presented for the k~c 
model calculations. These are U/u* in semi-log plot, U/U 0 , 
uv/u| , k/u^ and eR/uJ versus l-2r/D. 

The results for the full Reyno 1 ds-s tress closure using (5) 
for diffusion modelling are presented in Figures 1-3. Rotta’s 
model results are given in Figure la-lh, while Launder et al.’s 
model calculations without and with wall correction are given in 
Figures 2a-2h and 3a-3h, respectively. Results for full 
Reynolds-s tress closure using (6) for diffusion modelling are 
shown in Figures 4-6. The first eight figures give the 
calculated results of Rotta’s model and the other figures those 
of Launder et al.’s model without and with wall correction. 
Algebraic stress model calculations are presented in Figures 
7-10. The first sixteen figures give the results of the 
equilibrium turbulence calculations using Rotta’s and Launder 
et al . ’ s model. This is followed by the same calculations 
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assuming non-equilibrium turbulence. Finally, the k-e results 
are presented in Figures lla-lle. 

For the sake of clarity, a discussion of these results 
and their comparisons with measurements is presented in three 
different sections. These are: (1) effects of mean-strain 

modelling, (2) effects of diffusion modelling and (3) effects of 
wall correction on redistribution modelling. Finally, the k ~c 
equation results are presented in Section 5.4. 

5 . 1 Effects of Mean -Strain Modelling 

The equations (17)-(22) with boundary conditions (23) and 
(24) are solved assuming C lw = C 2w = 0. For each experiment, two 
calculations are carried out; one with a x = /S x - Y 1 ~ 0 and 

another with these constants as given in Table 1. Therefore, a 
comparison of these two calculations with the measured data would 
reveal the relative merits of mean-strain modelling. The results 
are presented in Figures 1 and 2. Other comparisons of the 
effects of mean-strain modelling, subject to different 
approximations for turbulent diffusion, are given in Figures 4 
and 5 for isotropic diffusion modelling, in Figures 7 and 8 for 
equilibrium algebraic stress closure and in Figures 9 and 10 for 
non-equilibrium algebraic stress closure. For ease of reference 
later on, the set of figures 1, 2; 4, 5; 7, 8 and 9, 10 shall be 

designated as set A, B, C and D, respectively. Therefore, figure 
sets A and B give the full Reyno 1 ds-s tress closure results, while 
sets C and D show the algebraic stress closure calculations. A 
comparison of the results within each set will indicate the 
relative merits of mean - strain modelling given a fixed diffusion 
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model for closure, while comparisons between different sets will 
elucidate the effects of diffusion modelling. 

In general, both redistribution models give reasonably good 
results for U, uv, k and c (Figures a-e of each set of Figures A, 
B, C, D). However, they fail to replicate correctly the 
behaviour of the normal stresses near the wall (Figures f-h of 
each set of Figures A, B, C, D). The calculated mean U is in 
good agreement with the measured mean U in the near wall region, 
but shows substantial discrepancy in the pipe core (Figures b of 
each set of figures), even though the logarithmic behaviour in 
this region is recovered (Figures a of each set). The measured U 
can be correlated by a logarithmic law-of- the-wal 1 , such that 

J- = «lnY + + p , (35) 

U * 

where a = 1 /k, k is the von Karman constant and p is 
parame t ical ly dependent on R D (Afzal and Yajnik 1973). For the 
experiments of Laufer (1954) and Schildknecht et al. (1979) the 
constants thus determined are listed in Table 2 together with 
the quantities U Q /u* . Likewise, these quantities can also be 
determined from the calculations. They are also listed in 
Table 2 for comparison. It can be seen that the measured slope 
of the log-law is not in agreement with the calculated slopes and 
that the calculated U Q /u* are always lower than the measured 
values. In view of this, it is very difficult to conclude which 
of the two redistribution models gives a better description of 
the mean flow. 
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In order to understand the discrepancy noted between the 
measured and calculated U, the mean U obtained by integrating 
(17) using the measured uv as input is also shown in Figures a of 
each set of figures for comparison with the model calculations. 
The corresponding a, p and U 0 /u* are listed in Table 2. It can 
be seen that the mean U thus determined is in excellent agreement 
with that calculated from Rotta’s model and the k -c closure. On 
the other hand, Launder et al.’s model consistently under- 
predicts the slope and U Q /u* but over-estimates the constant p . 
Both models, however, predict an increase in p as R D is decreased 
just as in the analysis of Afzal and Yajnik (1973). As for the 
behaviour of uv, k and e , the two model calculations are in 
excellent agreement with measurements, especially near the wall. 
This demonstrates that the modified dissipation model is valid 
and can account for the near wall behaviour very well. 

The performance of this dissipation model is quite 
independent of the redistribution model (Figure sets A and B). 
Therefore, based on the above comparison, it can be concluded 
that, as far as the mean flow, TTv , k and c behaviour's are 
concerned, the simple ret urn- t o- i so t ropy model of Rotta is just 
as promising as the more complete model of Launder et al. (1975). 
The performance of the Launder et al. (1975) model is found to be 
not as good when used in conjunction with algebraic stress 
closures (Figure sets C and D). In view of this, mean strain 
modelling is found to have a negative effect on the overall flow 
behaviour when turbulent diffusion is improperly modelled 
(compared figure sets A and B, C and D). A similar conclusion 
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has also been reached by Yao and So (1985) in their analysis of 
curved-pipe flow using rapid distortion theory. The present 
results together with Yao and So’s (1985) analysis, therefore, 
point to the modelling of the turbulent redistribution terms by 
return-to-isot ropy models for fully-developed pipe flows. 

Even though different levels of v 2 and w 2 are predicted by 
Launder et al.’s model, their comparison with measured data is 
slightly better than those obtained from Rotta's model (Figures f 
of each set). Essentially, both models under-estimate the rise 
of u 2 and greatly over-estimate the rise of v 2 and w 2 near the 
wall. They also fail to predict the isotropic behaviour of the 
turbulence field at the pipe center. As a result, the k 
distribution in the pipe core is over-predicted (Figures d of 
each set). A comparison of the stress ratios, u^^/k, clearly 
shows the inadequacy of the two models (Figures g and h of each 
set). Launder et al.’s model gives a better correlation with 
data for uv/k (Figures h of each set); however, it leads to a 
rather flat variation for u { u. /k (Figures g of each set). 

Besides, the limiting values of u i u i /k at the wall are not 
predicted correctly. In evaluating the limiting values, the 
measured data is fitted to the expansions 

u i u j 17 a io y N + b i j y N+1 + • • • 

near a wall. This allows the a d ^ and b i * to be determined and 
hence the values (u. u. /k) . For example, experimental values 

thus determined are: (u 2 /k) w ~ 2, (v 2 /k) w = 0 and (w 2 /k) w ^ 0 for 

Laufer’s data. The corresponding calculated values are 
(u 2 /k) w ^ .9, (v 2 /k) w = (w 2 /k) w ~ .55 for Rotta’s model and 
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(u 2 /k) w ^ .9, (v 2 /k) w ^ .5 and (w 2 /k) w .6 for Launder et al.’s 

model. Therefore, the rapid decrease of u 2 /k and steep rise of 
v 2 /k and w 2 /k near a wall are not predicted by the models at all 
(Figures g of each set). Rased on these calculations, it seems 
that the assumption, Uj u i /k are uniform across the flow, is 
inherent in these closure models. The effects of diffusion 
modelling on the behaviour of u A u^ /k near a wall will be further 
examined in Section 5.2. 

5 . 2 Effects of Diffusion Model ling 

In order to investigate the effects of diffusion modelling 
on turbulence closure, another set of calculations is performed 
with a vastly different diffusion model. This time an isotropic 
gradient diffusion model for the turbulent stresses is assumed as 
in (6). The diffusivity v t is taken to be given by (7) with the 
damping function f 2 near a wall included to account for wall 
proximity effects. Again, calculations are carried out for both 
the Rotta and Launder et al. models for the redistribution terms. 
The solution of (17)--(22) is performed assuming C lw and C 2w to be 
identically zero. 

The results are also plotted in Figures sets A-D for 
comparison with the previous calculations. Calculated values for 
a, p and U 0 /u* cire listed in Table 2. It can be seen that 
Rotta’ s model now gives rise to over-prediction of « and 
U 0 /u* and under-estimation of p. However, the parametric 
dependence of p on R r , is not correctly p r e d i c t e d . On the other' 
hand, Launder et al.’s model also gives an incorrect trend for p . 
Instead of predicting an increase for p as R D decreases, it gives 
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a that decreases slightly with R D . Also, the slope, cx, 
increases as R D decreases. In view of these incorrect trends, it 
can be concluded that the performance of Launder et al. *s model 
does not fair well with Rotta’s model when an isotropic gradient 
diffusion model is used to approximate turbulent diffusion. 

The near wall behaviours of U, uv , k and e are again well 
predicted. This shows that the modelled flow near a wall is 
essentially governed by the dissipation model and is only 
slightly dependent on the diffusion and pressure redistribution 
models. There are small differences in the calculations of 

u i u i/ u * and n i u./k (Figures f-h of each set). However, they are 
not significant enough to warrant a conclusion that one diffusion 
model is better than another. Essentially, the shortcomings 
noted in Section 5.1 for the diffusion model given by (5) are 
also true for (6). Therefore, the calculations indicate that 
once a gradient diffusion model is assumed, the results are only 
slightly dependent on the behaviour of the diffusivity. An 
isotropic model will give results that are quite similar to those 
obtained from a non-isotropic one. 

5 . 3 Effects of Wall Correction on Redistribution Modelling 

The effects of wall correction on redistribution modelling 
are assessed by solving (17)-(22) with the wall correction terms 
included. In the course of solving these equations, it is found 
that if the C lw and C 2w values suggested by Launder et al. (1975) 

are used convergent solutions to (17)-(22) are not possible. The 

2 

problem is traced to the coefficient of the terms (u i u i - ^k) 
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in (19)-(21). After normalization, the coefficient becomes 
[C 1 /Re - C lw (^ 3 + ) 3/2 /2 3/2 ^ 2 v], If C lw = . 125 is used, 

as suggested by Launder et al. , the coefficient becomes negative 

over a substantial portion of the pipe. Consequently, these 

2 

terms (u.u* - ^-k) change sign and the equations are not balanced. 

A similar behaviour is also observed in the production term in 
(20) because C 2w is too large. Subsequently, C lw and C 2w are 
slowly decreased until convergent solutions to (17)-(22) are 
obtained. The values of C lw and C 2w thus determined are shown in 
Table 1. Two sets of model calculations are performed and these 

are carried out with Launder et al.’s model for the 
redistribution term and (5) and (6) for the diffusion model. 

Actually, convergent solution is possible at some higher 
values of C lw and C 2w . However, the calculated results compare 
poorly with measurements and they are not shown. Calculations 
have been made with a series of values for C lw and C 2w and they 
lie in the range .025 < C lw < .0625 and .003 < C 2w < .0075. The 
largest values denote the upper limit for C lw and C 2w where 
convergent solution is possible. In general, the effects of wall 
corrections are to increase a, p and U 0 /u^, thus increasing the 
discrepancies noted between calculations and measurements of the 

mean flow. Furthermore, the peak value predicted for k increases 
as C lw and C 2w are increased. For example, the calculated peak 
value for k is more than 20% higher than the measured value when 
C lw = .0625 and C 2w = .0075 are used. Reduction of these 
constants to the values given in Table 1 gives the results shown 
in Figures 3 and 6. Even then, the peak value of k is . 10% 
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higher than measurements and the calculations obtained by 
neglecting wall corrections (compare Figures 2 and 3, 5 and 6). 

On the other hand, when C lw = .025 and C 2w = .003 are used in the 
governing equations, the calculated results are essentially 
indentical to those shown in Figures 2 and 5. If plotted, they 
practically overlap on top of the curves shown in Figures 2 and 
5. Even in the near wall region, little differences are noted. 

In view of these results, it can be said that if large values of 
C lw and C 2w are specified, wall corrections affect the 
calculations adversely. However, when small values of C lw and 
C 2w are used for the calculations, wall corrections have little 
effect on the results. 

5 . 4 Two-Equation Closure Results 

The results obtained by solving (14), (15) and (17) together 

with the appropriate boundary conditions are shown in Figure 11. 
They display characteristics very similar to those obtained from 
an equilibrium ASM closure using Rotta’s model (Figure 7) and 
from a full Reynolds-s tress closure using (5) for diffusion 
modelling and Rotta’s return-to-isot ropy model for pressure 
redistribution (Figure 1). The calculated a, p and U Q /u* for 
these three cases are very similar. However, both the k-c 
closure and the equilibrium ASM closure fail to predict the 
correct increase in p as R D is decreased (Table 1 and Afzal and 
Yajnik 1973). In view of this, the performance of the two- 
equation closure is not as good as that of the full Reynolds- 
stress closure. On the other hand, its prediction of fully- 
developed turbulent pipe flow properties are better than any ASM 
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closure models that used the Launder et al. model for pressure 
redistributions (compare Figure 11 with Figures 8 and 10 and the 
values of a, p listed in Table 2). Therefore, if only U, uv, k 
and e information are required in any pipe flow calculation, the 
k-e equation is a simple closure model to use and will provide 
reliable results. 
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6 . 


Application to Ful ly-Deve loped Rotating Pipe Flows 


One of the objectives of the present study is to develop a 
closure model for solid fuel ramjet combustor calculation. If 
the solid fuel ramjet combustor is used to power a projectile 
fired from a cannon, the whole combustor would spin at a very 
high rate. As a result of this spin, the flow inside the 
combustor would also be subject to the influence of a large 
circumferential velocity which has its maximum value at the 
combustor wall and decreases rapidly to zero at the combustor 
centerline. Therefore, the usual logarithmic law-of-the wall may 
not apply to the flow very near the wall and the high-Reynolds- 
number closure model may not be applicable to this kind of 
combustor flow calculation. With the development of the low- 
Reynolds-number closure models discussed in Section 5, they can 
then be applied to assess the effects of rotation on the 
calculated flow field correctly. This section presents the 
results of such an assessment. 

The effects of rotation can be best illustrated by 
considering a simple model problem where rotation appears as the 
only additional parameter in the flow field. Such a problem is 
given by the fully-developed flow through a circular pipe 
rotating at a constant speed. Since the above results indicate 
that the equilibrium ASM closure using Rotta’s model is just as 
good as a full Reynolds-s tress closure model, the following 
analysis is carried out with the equilibrium ASM closure model 
only. In Section 6.1, the governing equations for the rotating 
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pipe flow problem are specified. The results of this calculation 
are discussed in Section 6.2 


6 . 1 Governing Equations for Flows Through a Rotating Pipe 

Cylindrical coordinates are again used to analyze the flow. 
The pipe is assumed to rotate at a constant angular speed of Q, 
so that the circumferential velocity of the fluid at the pipe 
wall is W 0 = RQ. When the flow becomes fully-developed and 
axisymmetric , d/dO = 0, d/dx = 0 and V = 0. The resultant mean 
momentum equations reduce to: 


0 


1 dp 1 d f dlk 

p dx r dr^ dr' 



dF (ruv) 


(35) 


W n 1 d / ' 

r>J r 2 dr r " 





If the equilibrium ASM with Rotta’s return-t o-isot ropy model 
is used to determine uv and vw, then the component equations for 
uVu j can be obtained from (11) and are: 
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If it is 

further 

assumed 

that 

when fully-developed flow is 



established, a solid-body rotation exists in the fluid, then 

W = , (46) 


and it follows from (36) that vw = 0. With these simplifi- 
cations, (35) can be integrated to give (17) and (37) and (38) 
reduce to (14) and (15) respectively* The solution of (40) to 
(45) then gives 


V 2 = W 2 = _ k , 


-rr C, - 1 2, 2k — dU 


u 


2 - 


3 k “ cTT uv d? ’ 


— 2 k 
” UW = C77 UV R^ * 


— C. - 1 2 k 2 

UV — 1 

C 2 3 c 


1 +^ k2 -A) 2 

C^cs'-R ‘ 


dU 

dr 


(47) 

(48) 

(49) 

(50) 
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Therefore, W Q influences uv according to (50) and uv in turn 
affects U , u 2 , k and c . The problem of ful ly-devel oped turbulent 
flow in a rotating pipe is described by equations (14), (15), 

(17) and (47) - (50). These equations are solved by the same 
technique discussed in Section 4. The boundary conditions are 
given by (39) rather than by (23) and (24). 

6 . 2 R esul ts 

Since there are no measurements available for comparison 
with the present calculations, only parametric studies are 
carried out. In order to evaluate the effects of rotation on the 
turbulence field, the rotation calculations should be carried out 
with a known condition for the non-rotating case. Therefore, the 
Laufer and Schildknecht et al experiments are selected as the 
known non-rotating case and parametric studies on the effects of 
rotation are carried out with these cases as the base. Three 
different calculations are performed. These are W Q = .105U Q , 

.21U 0 and .42U 0 , and the corresponding Q are 24 RPM, 48 RPM and 
96 RPM, respectively. The results are shown in Figure 7 for 
comparison with the zero rotation case. 

In general, rotation has a great influence on the flow even 
in the fully-developed state. The effects of rotation on the 
mean flow are clearly evident in the pipe core (Figures 7a and 
7b) and it tends to decrease the extent of the log-law region as 
Q increases. Also, the turbulent kinetic energy in the pipe core 
is increased by flow rotation because mixing is being promoted 
due to the action of the centrifugal forces. At very high 
rotation, i.e. Q - 96 RPM, k remains fairly constant in the pipe 
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core after the maximum value is reached near the pipe wall 
(Figure 7d) . Similar trends are also noted for other turbulence 
properties (Figures 7e-g). However, since uv remains unaffected 
by Q in a fully-developed flow (Figure 7c), uv/k decreases with 
rotation (Figure 7h) and the region where uv/k is constant 
disappears once there is rotation in the pipe. This shows that 
the assumption, uv/k = constant over a substantial portion of the 
pipe, normally invoked by turbulence modellers for simple 
turbulent flows is not valid for rotating turbulent flows. 
Finally, the shear stress uw is not small, and dependent on Q, 
can even be larger than uv (Figure 7i). The maximum of uw does 
not occur very near the pipe wall as in the case of uv . It 
occurs about half-way between the wall and the pipe center. 


41 



7. 


Conclusions 


A low-Reynolds-number turbulence closure for the full set of 
Reynolds-s tress equations is formulated. The formulation is 
based on a gradient diffusion model for turbulent diffusion, 
conventional high-Reynolds-number model for pressure 
redistribution and a modified dissipation model that accounts for 
viscosity effects near a rigid wall. Validation of the closure 
model is carried out with fully-developed turbulent pipe flows at 
two different Reynolds number. In general, the models give good 
results for U, uv , k and r, but fails to reproduce the behaviour 
of the normal turbulent stresses. The failings of the model are 
even more evident when the structure parameters, u.Uj/k, are 
compared and are especially noticeable in the near wall region. 
These discrepancies cannot be erased by modifying the 
redistribution model to account for the reduced turbulence energy 
transfer from the streamwise direction to that normal to a wall 
as suggested by Launder et al. (1975). Neither can the 
correlations between pre diction and measurement s be improved by 
the inclusion of mean-strain terms in the modelling of the 
pressure redistribution terms. Furthermore, it is found that 
once the gradient diffusion assumption is invoked, the calculated 
results are only slightly dependent on the diffusivity 
assumption. An isotropic diffusivity model will give results 
that are quite similar to those obtained by using a non- isotropic 
one . 

All full -Reynol ds-s tress closure models examined show the 
same shortcomings when applied to calculate ful ly -developed 
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turbulent pipe flows. They all fail to predict the steep rise of 
the normal stresses near the pipe wall and the isotropic 
behaviour of the turbulence field at the pipe center. However, 
the Launder et al model for pressure redistribution seems to 
provide better agreement with data concerning the prediction of 
the mean flow. In view of these results, it can be concluded 
that a closure model based on the Launder et al. model for the 
pressure redistribution and a non-isotropic gradient diffusion 
model gives the best overall results for fully-developed 
turbulent pipe flow calculations. 

The same conclusion cannot be reached when algebraic stress 
closures are considered, however. Here, Launder et al.*s model 
gives results that are less appealling than those obtained with 
Rotta’s model. As before, the manner in which turbulent 
diffusion is modelled has little effect on this overall 
conclusion; that is, it is true for equilibrium ASM as well as 
for non-equilibrium ASM. The algebraic stress closures give 
results that are closely similar to those obtained from full 
Reynolds-stress closures. Only minor differences appear in the 
predictions of the stress ratios, u.u^/k. Otherwise the 
performance of the ASM closures is just as good. The same can 
also be said of the two-equation closure. Even then, the best 
overall prediction of the logarithmic law-of- the-wal 1 is provided 
by the full Reynolds stress model using a non-isotropic diffusion 
model and the Launder et al. model for pressure redistribution. 

Finally, a parametric study of rotation effects on fully- 
developed turbulent pipe flows reveals that increasing rotation 
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decreases the extent of the log-law region. Also, rotation tends 
to increase the overall level of turbulent kinetic energy in the 
pipe core, as well as other turbulent properties. 
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Table 1 

Model constants and damping 

functions 
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Table 2. A Comparison of the Calculated and Measured 
Constants in the Logarithmic Law-of-the-wal 1 


■ 

Laufer (1954) 

Schildknecht 
et al. (1979) 

• 

a 

<*) 

P 


a 

(*) 

P 

U 0 / U * 

Measured U 

2.50 

(0.40) 

5.20 

23.76 

2.50 

(.40) 

5.45 

22.25 

U from measured 
uv 

2.60 

(0.385) 

5.00 

22.29 

2.53 

(0.395) 

5.90 

21.23 


Rotta* s 
model 

2.65 

(0.377) 

5.00 

22.69 

2.66 

(0.378) 

5.31 

21.08 

Non-isotropic 

diffusion 

model 

Launder 
et al . * s 
model 

2.46 

(0.406) 

5.69 

22.10 

2.48 

(0.403) 

6.00 

20.79 

Launder 
et al . 

+ wall 
correction 

2.62 

(0.382) 

5.53 

23.04 

2.82 

(0.355) 

4.97 

21.65 

Isotropic 

diffusion 

model 

Rotta’ s 
model 

2.62 

(0.382) 

5.00 

22.60 

2.69 

(0.372) 

5.00 

20.79 

Launder 
et al’s 
model 

2.35 

(0.426) 

6. 10 

21.70 

2.45 

(0.408) 

5.95 

20.43 


Launder 
et al * s 
+ wall 
correction 

2.30 

(0.435) 

6.97 

22.40 

2.42 

(0.414) 

6.72 

21.06 

Non- 

equilibrium 

turbulence 

Rotta* s 
model 

2.56 

(0.391) 

6 

23.21 

2.65 

(0.377) 

5.81 

21.55 

Launder 
et al . * s 
model 

2.40 

(0.417) 

6.3 

22.04 

2.42 

(0.413) 

6.3 

20.45 


Rotta’ s 
model 

2.67 

(0.375) 

5.06 

22.77 

2.64 

(0.379) 

5.36 

21.15 

Equilibrium 

turbulence 

Launder 
et al . * s 
model 

2.47 

(0.405) 

4.63 

21.27 

2.51 

(0.398) 

4.78 

19.78 

k -e eq. 
closure 

2.62 

(0.382) 

5.06 

22.76 

2.64 

(0.379) 

5.36 

21.15 
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Figure la. Comparison of Laufer's (1954) and 
Schildknecht et al.'s (1979) measurements with model 
calculations. Model is based on full Reynolds -stress 
closure with (5) and (10) used for diffusion and 
dissipation modeling. Pressure redistribution is 
provided by Rotta's (1951) model. 
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Figure 2a. Cornparisoii -of. Laufer's.. (1954) and Schildknecht 
et al.'s (1979) measurements with model calculations. 
Model is based on full Reynolds -stress closure with (5) 
and (10) used for diffusion and dissipation modelling. 
Pressure redistribution is provided by Launder et al.'s 
(1975) model given in (3). 


Ln Y + 












k / u * / u * 





R/u*/u*/u* 


o 


o_ J 

10 7 / 


~2r- — __ 


0.00 


20.00 


40.00 

Y + 


60.00 


Q LAUFERI 1954 ) 

&■ SCHILDKNECHT et a l . { 1979) 

CALCULAT! ONC RaD-500CC ! 

• CALCl)LAT!0N(ReP”4!750) 


Figure 2e 














LAUFER! 1954 ) 

SCH I LDKNECHT at al .( 1979 ) 
CALCULATION! RaD-50000 ) 
CALCULATION! RaD-21750 ) 







U/u* 


EXP. DATA - LAUFERC I9S4 ) 

U FROM MEASURED -uv - LAUFERM9S4) 

EXP. DATA - SCH I LDKNECHT at q|.M979) 

U FROM MEASURED -uv - SCH I LDKNECHT at al.(1979) 
CALCULAT 1 0 N( RcD-50000 S 
CALCULATION! ReD-21750 ) 


A Q 


A A 

A 
A — 



Figure 3a. Comparison of Laufer's (1954) and 
Schildknecht et al.'s (1979) measurements with model 
calculations. Model is based on full Reynolds -stress 
closure with (5) and (10) used for diffusion and 
dissipation modelling. Pressure redistribution with 
wall correction is provided by (4) . 
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Figure 4a. Comparison o£ Lau£er T s (1954) and Schildknecht et al. 
(1979) measurements with model calculations. Model is based on 
full Reynolds-stress closure with (6) and (10) used for diffusion 
and dissipation modelling. Pressure redistribution is provided 
by Rotta's (1957) model. 
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Figure 5a. Comparison of Laufer's (1954) and 
Schildknecht et al.'s (1979) measurements with model 
calculations. Model is based on full Reynolds -stress 
closure with (6) and (10) used for diffusion and 
dissipation modelling. Pressure redistribution is 
provided by Launder et al.'s (1975) model given in (3) 
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Figure 6a. Comparison of Laufer's (1954) and Schildknecht et 
al.'s (1979) measurements with model calculations. Model is 
based on full Reynolds -stress closure with (6) and (10) used 
for diffusion and dissipatation modelling. Pressure redistri 
bution with wall correction is provided by (4) . 
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Figure 7a. Comparison of Laufer's (1954) and 

Schildknecht et al. 's (1979) measurements with model 
calculations including rotation. Model is based on 
ASM closure provided by (11), (12), (14) and (15). 
Pressure redistribution model is given by Rotta 
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Figure 8a. Comparison o£ Laufer’s (1954) and 

Schildknecht et al. 's (1979) measurements with model 
calculations. Model is based on ASM closure provided 
by (11) , (12) , (14) and (15) . Pressure redistribu- 
tion model (3) is given by Launder et al. (1975). 
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Figure 9a. Comparison of Laufer's (1954) and Schildknecht 
et al.'s (1979) measurements with model calculations. 

Model is based on ASM closure provided by (12) -(15), 
Pressure redistribution model is provided by Rotta (1951) 
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Figure 10a. Comparison o£ Laufer's (1954) and Schildknecht 
et al.'s (1979) measurements with model calculations. 
Model is based on ASM closure provided by (12) -(15). 
Pressure redistribution model (3) is given by Launder 
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Figure 10b. 
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Figure 11a. Comparison o£ Laufer's (1954) and Schildknecht 
et al.'s (1979) measurements with model calculations. 
Model is based on two -equation closure provided by 
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